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These lectures provide an overview of Quantum Chromodynamics (QCD), the 
SU{3)c gauge theory of the strong interactions. The running of the strong coupling 
and the associated property of Asymptotic Freedom are analyzed. Some selected 
experimental tests and the present knowledge of cts and the quark masses are 
briefly summarized. A short description of the QCD flavour symmetries and the 
dynamical breaking of chiral si/mm«i*j*jp| also given. A morq-iAstailed discussion 
can be found in standard textbooks El'B'Qo and recent reviews Q'Q. 



1 Quarks and Colour 

A fast look into the Particle Data Tables reveals the richness and variety 
of the hadronic spectrum. The large number of known mesonic and baryonic 
states clearly signals the existence of a deeper level of elementary constituents 
of matter: quarksu. In fact, the messy hadronic world can be easily understood 
in terms of a few constituent spin-^ quark flavours: 



g = +i 


u c t 




d s b 



Assuming that mesons are M = qq states, while baryons have three quark 
constituents, B = qqq, one can nicely classify the entire hadronic spectrum. 
There is a one-to-one correspondence between the observed hadrons and the 
states predicted by this simple classification. Thus, the Quark Model appears 
to be a very useful Periodic Table of Hadrons. 

However, the quark picture faces a problem concerning the Fermi-Dirac 
statistics of the constituents. Since the fundamental state of a composite sys- 
tem is expected to have L — 0, the A++ baryon (J = |) with J3 = +| 
corresponds to u^u^u^ , with the three quark spins aligned into the same di- 
rection (s3 = and all relative angular momenta equal to zero. The wave 
function is symmetric and, therefore, the A''"'' state obeys the wrong statis- 
tics. The problem can be solved assumingd the existence of a new quantum 
number, colour, such that each species of quark may have Nc — 3 different 
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Figure 1: Two-jet event from the hadronic decay of a Z boson (DELPHI). 



colours: g", a — 1,2,3 (red, yellow, violet). Then, one can reinterpret this 
state as ~ e"^'^ \ul^Upul^) (notice that at least 3 colours are needed 

for making an antisymmetric state). In this picture, baryons and mesons are 
described by the colour-singlet combinations 

S = _L ,"/37 Iq^q^q^^ , ^'^ = ^ ^'^^ ''^"^^^ ' ^^'^^ 

In order to avoid the existence of non-observed extra states with non-zero 
colour, one needs to further postulate that all asymptotic states are colourless, 
i.e. singlets under rotations in colour space. This assumption is known as the 
confinement hypothesis, because it implies the non-observability of free quarks: 
since quarks carry colour they are confined within colour-singlet bound states. 

The quark picture is not only a nice mathematical scheme to classify the 
hadronic world. We have strong experimental evidence of the existence of 
quarks. Fig. |l| shows a typical Z — > hadrons event. Although there are 
many hadrons in the final state, they appear to be coUimated in 2 jets of 
particles, as expected from a two-body decay Z — s- qq, where the qq pair has 
later hadronized. 

1.1 Evidence of Colour 

A direct test of the colour quantum number can be obtained from the ratio 



Figure 2: e+e~ — ► hadrons. Figure 3: t decay. 

The hadronic production occurs through e+e^ j*,Z* qq ^ hadrons. 
Since quarks are assumed to be confined, the probability to hadronize is just 
one; therefore, summing over all possible quarks in the final state, we can 
estimate the inclusive cross-section into hadrons. At energies well below the Z 
peak, the cross-section is dominated by the 7-exchange amplitude; the ratio 
Re+e- is then given by the sum of the quark electric charges squared: 

Nf { = 2, {Nf = 3 : u,d,s) 



^cY.Q) ^ \f^c^f. iNf^4:u,d,s,c) . (1.3) 

lf^c = T' (^/ = 5 : u,d,s,c,b) 



3 
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The measured ratio is shown in Fig. Although the simple formula ( [l.3[ ) 
cannot explain the complicated structure around the different quark thresholds, 
it gives the right average value of the cross-section (away from thresholds), 
provided that Nc is taken to be three. The agreement is better at larger 
energies. Notice that strong interactions have not been taken into account; 
only the confinement hypothesis has been used. 

Additional evidence for Nc = 3 is obtained from the hadronic decay of the 
T lepton, which proceeds through the VF-emission diagram shown in Fig. ^ 
Since the W coupling to the charged current is of universal strength, there 
are (2 -I- Nc) equal contributions (if final masses and strong interactions are 
neglected) to the r decay width. Two of them correspond to the leptonic decay 
modes — > Vre^De and —t v-rfJ- i^tn while the other Nc are associated 
with the possible colours of the quark-antiquark pair in the — + Vrdeu decay 
mode {dg = cos Bed + sinOcs). Hence, the branching ratios for the different 
channels are expected to be approximately: 

Br^i = Br(T- ^ u^l-Di) « = i = 20% , (1.4) 

2 + i\c o 

T(t~ Vt ^ hadrons) 

^ - r^^«A^c = 3, 1.5 

1 (r v^e v^) 
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Figure 4: Measurements [J of 



The two continuous curves are QCD fits. 



which should be compared with the experimental averagesi 



0: 



B^^e = (17.791 ±0.054)%, 



Br^t, = (17.333 ± 0.054)% , (1.6) 



Rt — (1 — 



B^^f,)/Br^^ = 3.647 ± 0.014 . 



(1.7) 



The agreement is fairly good. Taking Nc — 3, the naive predictions only 
deviate from the measured values by about 20%. Many other observables, 
such as the partial widths of the Z and bosons, can be analyzed in a 
similar way to conclude that Nc = 3. 

A particularly strong test is obtained from the tt*' — s- 77 decay, which 
occurs through the triangular quark loops in Fig. ||. The crossed vertex denotes 
the axial current = {u^f^j^u — d-y^j^d). One gets: 



(1.8) 



where the tt" coupling to A^, /^r — 92 A MeV, is known from the tt h 
decay rate (assuming isospin symmetry). The agreement with the measured 
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Figure 5: Triangular quark loops generating the decay tt' 
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valueLO, r = 7.7 ± 0.6 eV, is remarkable. With Nc = 1, the prediction would 
have failed by a factor of 9. The nice thing about this decay is that it is 
associated with a quantum anomaly: a global flavour symmetry which is bro- 
ken by quantum effects (the triangular loops). One can proci-that the decay 



Anomalies provide another compelling theoretical reason to adopt Nc = 3. 
The gauge symmetries of the Standard Model of electroweak interactions have 
also anomalies associated with triangular fermion loops (diagrams of the type 
shown in Fig. H, but with arbitrary gauge bosons — , Z, 7 — in the external 
legs and Standard Model fermions in the internal lines). These gauge anomalies 
are deathly because they destroy the renormalizability of the theory. Fortu- 
nately, the sum of all possible triangular loops cancels if Nc — 3. Thus, with 
three colours, anomalies are absent and the Standard Model is well defined. 

1.2 Asymptotic Freedom 



The proton is an extended object with a size of the order of 1 fm. Its 
structure can be investigated through the electromagnetic scattering e^p — > 
e~p. At very low energies, the photon probe is unable to get information on 
the proton structure and the proton behaves as a pointlikc particle. At higher 
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Figure 6: Inelastic e p — > e X scattering. 
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energies, the photon is sensitive to shorter distances; the proton finite size 
gives then rise to form factors, which suppress the elastic cross-section at large 
values of = —q^ = AEE' sin.^ |, where q'^ = (fee — k'^)'^ is the momentum 
transfer through the intermediate photon propagator, E and E' the energies 
of the incident and scattered electrons in the proton rest frame, and 9 the 
scattering angle. 

One can try to further resolve the proton structure, by increasing the 
incident energy. The inelastic scattering e~p — *■ e~ X becomes then the dom- 
inant process. Making an inclusive sum over all hadrons produced, one has 
an additional kinematical variable corresponding to the final hadronic mass, 
W'^ = Px- The scattering is usually described in terms of and 

Mp 2Mp ' ^ ' 

where and Mp are the proton cuadrimomentum and mass; v is the en- 
ergy transfer in the proton rest frame. In the one-photon approximation, the 
unpolarized differential cross-section is given by 

da Tra^cos^l ( , 9 , .9^ 96*1 

The proton structure is then characterized by two measurable structure func- 
tions. For a pointlike proton, the elastic scattering corresponds to 



^^iQ'^^)-^/[^-it;)^ (1.11) 



4 

At low Q^, the experimental data reveals prominent resonances; but this 
resonance structure quickly dies out as increases. A much softer but size- 
able continuum contribution persists at large Q^, suggesting the existence of 
pointlike objects inside the proton. 

To get an idea of the possible behaviour of the structure functions, one 
can make a very rough model of the proton, assuming that it consist of some 
number of pointlike spin-i constituents (the so-called partons), each one car- 
rying a given fraction £j of the proton momenta, i.e. Pi^ = ^iP^. That means 
that we are neglecting the transverse parton momenta, and ~ £,Mp. The 



These approximations can be made more precise going to the infinite momentum frame of 
the proton, where the transverse motion is negligible compared with the large longitudinal 
boost of the partons. 
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interaction of the photon probe with the parton i generates a contribution to 
the structure functions given by: 



W^'>iQ\,.) = eiS[u-^^ef-6{C.-x), (1.12) 



where et is the parton electric charge and 



2MpV Q^ + W^- Af2 



(1.13) 



Thus, the parton structure functions only depend on the ratio x, which, more- 
over, fixes the momentum fractions ^i. We can go further, and assume that in 
the limit — > cxd, cx3, but keeping x fixed, the proton structure functions 
can be estimated from an incoherent sum of the parton ones (neglecting any 
strong interactions among the partons). Denoting fi{S.i) the probability that 
the parton i has momentum fraction ^i, one then has: 

, p ■ p 



W2{Q\ = L /*(^') w^\q^' ^) ^ ^ E /^(^) ^ I ^2(2:) . (1.14) 



This simple parton description implies then the so-called BjorkenEl scaling: the 
proton structure functions only depend pn the kinematical variable x. More- 
over, one gets the Callan-Gross relationEj 

F2{x) ^ 2xFi{x) , (1.15) 

which is a consequence of our assumption of spin- ^ partons. It is easy to check 
that spin-0 partons would have lead to Fi{x) = 0. 

The measured values of i^W2{Q^, v) are shown in Fig. ^ as function of x, 
for many different values of between 2 and 18 GeV^; the concentration of 
data points along a curve indicates that Bjorken scaling is correct, to a quite 
good approximation. Fig. ^ shows that the Callan-Gross relation is also well 
satisfied by the data, supporting the spin-i assignment for the partons. 

The surprising thing of these successful predictions is that we have assumed 
the existence of free independent pointlike partons inside the proton, in spite of 
the fact that quarks are supposed to be confined by very strong colour forces. 
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Figure 7: Experimental data on i'W2 
as function of x, for different values of 
(taken from Ref. 1). 




Figure 8: The ratio 2afft/^2 versus x, 
for different values llf (1.5 GeV^ < 
< 16 GeV^) (taken from Ref. 1). 



Bjorken scaling suggests that the strong interactions must have the property 
of asymptotic freedom: they should become weaker at short distances, so that 
quarks behave as free particles for — > oo. This also agrees with the empirical 
observation in Fig. ^, that the free-quark description of the ratio i?e+e- works 
better at higher energies. 

Thus, the interaction between a qq pair looks like some kind of rubber 
band. If we try to separate the quark form the antiquark the force joining 
them increases. At some point, the energy on the elastic band is larger than 
2m,qi , so that it becomes energetically favourable to create an additional q'q' 
pair; then the band breaks down into two mesonic systems, qq' and q'q, each 
one with its corresponding half-band joining the quark pair. Increasing more 
and more the energy, we can only produce more and more mesons, but quarks 
remain always confined within colour-singlet bound states. Conversely, if one 
tries to approximate two quark constituents into a very short-distance region, 
the elastic band loses the energy and becomes very soft; quarks behave then 
as free particles. 

1.3 Why SU{i)? 

Flavour-changing transitions have a much weaker strength than processes me- 
diated by the strong force. The quark flavour is associated with the electroweak 
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interactions, while strong forces appear to be flavour-conserving and flavour- 
independent. On the other side, the carriers of the electroweak interaction (7, 
Z, W^) do not couple to the quark colour. Thus, it seems natural to take 
colour as the charge associated, with the strong forces and try to build a quan- 
tum field theory based on it t3. The empirical evidence described so far puts 
a series of requirements that the fundamental theory of colour interactions 
should satisfy: 

1. Exact colour symmetry Gc (hadrons do not show colour multiplicity). 

2. Nc — 3. Thus, quarks belong to the triplet representation 3 of Gc- 

3. Quarks and antiquarks are different states. Therefore, 3* 7^ 3. 

4. Confinement hypothesis: hadronic states are colour singlets. 

5. Asymptotic freedom. 

Among all compact simple Lie groups there are only four having irreducible 
representations of dimension 3; three of them are isomorphic to each other. 
Thus, we have only two choices: SU{3) or SO{3) ~ SU{2) ~ Sp{l). Since 
the triplet representation of 5*0(3) is real, only the symmetry group SU{3) 
survives the conditions 1, 2 and 3. The well-known SU{3) decomposition of 
the products of 3 and 3* representations, 

qq: 3(g)3* = 1©8, 

qqq: 3 (g) 3 (g) 3 = 1 ® 8 8 10 , 

qq : 3(g3 = 3* ©6, 

qqqq : 3(g)3(g)3(g)3 = 3e3e3e6*e6*©15el5@15@15', 

guarantees that there are colour-singlet configurations corresponding to meson 
(qq) and baryon {qqq) states, as required by the confinement hypothesis. Other 
exotic combinations such as diquarks {qq, qq) or four-quark states {qqqq, qqqq) 
do not satisfy this requirement. 

Clearly, the theory of colour interactions should be based on SU{3)c- It 
remains to be seen whether such a theory is able to explain confinement and 
asymptotic freedom as natural dynamical consequences of the colour forces. 

2 Gauge Symmetry: QED 

Let us consider the Lagrangian describing a free Dirac fermion: 

£0 = i'^(x)7^5p*(x) - m¥(x)*(a;). (2.1) 
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Cq is invariant under global U{\) transformations 



^-(2;) ^ W{x) = exp{iQ6l}^'(a;), (2.2) 

where Q9 is an arbitrary real constant. The phase of ^(x) is then a pure 
convention-dependent quantity without physical meaning. However, the free 
Lagrangian is no-longer invariant if one allows the phase transformation to 
depend on the space-time coordinate, i.e. under local phase redefinitions 9 — 
9{x), because 

^ exp{iQ0} {d^ + iQd^e) -^{x) . (2.3) 

Thus, once an observer situated at the reference point a;o has adopted a given 
phase convention, the same convention must be taken at all space-time points. 
This looks very unnatural. 

The "Gauge Principle" is the requirement that the U{1) phase invariance 
should hold locally. This is only possible if one adds some additional piece 
to the Lagrangian, transforming in such a way as to cancel the d^O term in 



Eq. (2.3). The needed modification is completely fixed by the transformation 
( ^■3[ ): one introduces a new spin-1 (since d^O has a Lorentz index) field A^,(x), 
transforming as 

A^{x) ^ A'^{x) ^ A^{x) + -^d,e, (2.4) 
and defines the covariant derivative 



D^-^{x) = [d^ - ieQA^{x)] *(x) , (2.5) 
which has the required property of transforming like the field itself: 

D^-^{x) ^ {D^,-^)'{x) = e^^{iQ6}D^-^{x). (2.6) 
The Lagrangian 

C = i *(2;)7^L>^*(a;) - m¥(a;)*(a;) = £0 + eQyl,,(x) *(a;)7^*(x) (2.7) 

is then invariant under local U{1) transformations. 

The gauge principle has generated an interaction between the Dirac spinor 
and the gauge field , which is nothing else than the familiar vertex of Quan- 
tum Electrodynamics (QED). Note that the corresponding electromagnetic 
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charge eQ is completely arbitrary. If one wants yl^ to be a true propagating 
field, one needs to add a gauge-invariant kinetic term 

^Kin = —^F^^F^^, (2.8) 

where = d^^A^ — dyA^ is the usual electromagnetic field strength. A 
possible mass term for the gauge field, ^^m^A^^A^, is forbidden because it 
would violate gauge invariance; therefore, the photon field is predicted to be 



massless. The total Lagrangian in (2/7) and (2^) gives rise to the well-known 
Maxwell equations. 

From a simple gauge-symmetry requirement, we have deduced the right 
QED Lagrangian, which leads to a very successful quantum field theory. Re- 
member that QED predictions have been tested to a very high accuracy. 

3 The QCD Lagrangian 

Let us denote q'j a quark field of colour a and flavour /. To simplify the 
equations, let us adopt a vector notation in colour space: = {qj , qj , q'j) . 
The free Lagrangian 

^o^^^f {h'^d^, - rrif) qf (3.1) 
/ 

is invariant under arbitrary global SU {3)c transformations in colour space, 

qf — > (qf)' ^U°'pq'^^ , UW = U^U^1, det iJ = 1 . (3.2) 

The SU{3)c matrices can be written in the form 

[/ = exp |-i5sy6'a| , (3.3) 

where A" (a = 1, 2, . . . , 8) denote the generators of the fundamental represen- 
tation of the SU(3)c algebra, and 0a are arbitrary parameters. The matrices 
A'' are traceless and satisfy the commutation relations 

[A°,A^] = 2^r^'=A^ (3.4) 

with /"''^ the SU{3)c structure constants, which are real and totally antisym- 
metric. Some useful properties of SU{3) matrices are collected in Appendix A. 

As in the QED case, we can now require the Lagrangian to be also invariant 
under local SU(3)c transformations, 0a = 0a{x). To satisfy this requirement. 
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we need to change the quark derivatives by covariant objects. Since we have 
now 8 independent gauge parameters, 8 different gauge bosons G^i^{x), the 
so-called gluons, are needed: 



qj ^[d^-ig,G^{x)]qs. (3.5) 



Notice that we have introduced the compact matrix notation 



[G^[x)U^[^\ G>i{x). (3.6) 

' a(3 



We want D^^qf to transform in exactly the same way as the colour-vector g/; 
this fixes the transformation properties of the gauge fields: 

— > {Df^y ^UDt'U'' ; G'' — > (G^)' :^UGf'U'< - — (d'^U) . 

gs 

Under an infinitesimal SU{3)c transformation, 



Qf {q'}y = qf~i9s ) ^^-1f 

G^a (G^y ^G^~ d>^{6ea) + gsP'^Se, G^, . (3.7) 

The gauge transformation of the gluon fields is more complicated that the one 
obtained in QED for the photon. The non-commutativity of the SU{3)c ma- 
trices gives rise to an additional term involving the gluon fields themselves. For 
constant 69a, the transformation rule for the gauge fields is expressed in terms 
of the structure constants f thus, the gluon fields belong to the adjoint 
representation of the colour group (see Appendix A) . Note also that there is a 
unique SU{3)c coupling gs- In QED it was possible to assign arbitrary elec- 
tromagnetic charges to the different fermions. Since the commutation relation 
( |3.4[ ) is non-linear, this freedom does not exist for SU{S)c- 

To build a gauge-invariant kinetic term for the gluon fields, we introduce 
the corresponding field strengths: 

G'^'ix) = — [I?^, D''] = d^G" - d^Gf" - igs [G^, G"] = — G^^x) , 

gs ' 2 

Gr(x) = d^^G: ~ d^G'i + gsr'-GtG: . (3.8) 

Under a gauge transformation, 

G^"" — > (G^'^)' = UG^"" , (3.9) 



12 



and the colour trace Tr(G^''G^i/) = IG^'^G^j, remains invariant. 

Taking the proper normaUzation for the gluon kinetic term, we finally have 
the SU{3)c invariant QCD Lagrangian: 

^QCD ^ - J Grc;. + E 97 (i^D^ - ^f) <lf ■ (3-10) 

/ 

It is worth while to decompose the Lagrangian into its different pieces: 

f 

+ .9.GJ^5]g77^f^) (3.11) 

f \ / af} 

2 

rabc /nAi/^i^ Pi^ r^c 9s fabcr /-t^ /^v /-id /-te 

The first line contains the correct kinetic terms for the different fields, which 
give rise to the corresponding propagators. The colour interaction between 
quarks and gluons is given by the second line; it involves the SU{3)c matri- 
ces A°. Finally, owing to the non-abelian character of the colour group, the 
G^^^G"^ term generates the cubic and quartic gluon self-interactions shown in 
the last line; the strength of these interactions is given by the same coupling 
gs which appears in the fermionic piece of the Lagrangian. 



In spite of the rich physics contained in it, the Lagrangian (3.10) looks very 
simple because of its colour symmetry properties. All interactions are given 
in terms of a single universal coupling g^, which is called the strong coupling 
constant. The existence of self-interactions among the gauge fields is a new 
feature that was not present in QED; it seems then reasonable to expect that 
these gauge self-interactions could explain properties like asymptotic freedom 
and confinement, which do not appear in QED. 

Without any detailed calculation, one can already extract qualitative phys- 
ical consequences from >Cqcd- Quarks can emit gluons. At lowest order in g^, 
the dominant process will be the emission of a single gauge boson. Thus, the 
hadronic decay of the Z should result in some Z — > qqG events, in addition to 
the dominant Z qq decays discussed before. Fig. ^ clearly shows that 3-jet 
events, with the required kinematics, indeed appear in the LEP data. Similar 
events show up in e^e~ annihilation into hadrons, away from the Z peak. 

In order to properly quantize the QCD Lagrangian, one needs to add to 
'Cqcd the so-called Gauge-fixing and Faddeev-Popov terms. Since this is a 
rather technical issue, its discussion is relegated to Appendix B. 
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Figure 9: Three-jet event from the hadronic decay of a Z boson (DELPHI). 

4 Quantum Loops 

The QCD Lagrangian is rather economic in the sense that it involves a single 
coupling gs- Thus, all strong-interacting phenomena should be described in 
terms of just one parameter. At lowest order in gg (tree level), it is straight- 
forward to compute scattering amplitudes involving quarks and gluons. Un- 
fortunately, this exercise by itself does not help very much to understand the 
physical hadronic world. First, we see hadrons instead of quarks and gluons. 
Second, we have learnt from experiment that the strength of the strong forces 
changes with the energy scale: the interaction is very strong (confining) at low 
energies, but quarks behave as nearly free particles at high energies. Obviously, 
we cannot understand both energy regimes with a single constant gs, which 
is the same everywhere. Moreover, if we neglect the quark masses, the QCD 
Lagrangian does not contain any energy scale; thus, there is no way to decide 
when the energy of a given process is large or small, because we do not have 
any reference scale to compare with. 

If QCD is the right theory of the strong interactions, it should provide 
some dynamical scale through quantum effects. 

4-1 Regularization of Loop Integrals 

The computation of perturbative corrections to the tree-level results involves 
divergent loop integrals. It is then necessary to find a way of getting finite 
results with physical meaning from a priori meaningless divergent quantities. 
Let us consider the self-energy gluon loop in Fig. The correspond- 
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Figure 10: Gluon self-energy diagram. 

ing contribution in momentum space can be easily obtained, using standard 
Feynman rules techniques: 

The result is proportional to 5^, because there are two qqG vertices, and there 
is a SU{3)c factor, Tp = \, coming from the colour trace iTr(A°A'') = d'^^Tp. 

The problem appears in the momentum integration, which is clearly di- 
vergent J d'^k{l/k^) = Qo]. We can define the momentum integral in many 
different (and arbitrary) ways. For instance, we could introduce a cut-off M, 
such that only momentum scales smaller than M are integrated; obviously, 
the resulting integral would be an increasing function of M. Instead, it has 
become conventional to define the loop integrals through dimensianal regular- 
ization: the calculation is performed in D — 4 + 2e dimensions EZI. For e ^ 
the resulting integral is well defined: 

r d^k k'^ik-qf _ -i 
J (2^ k'^{k~qf ^ 6(4^ 

The ultraviolet divergence of the loop appears at e = 0, through the pole of 
the Gamma function, 

r(-e) = -i - 7s + 0{t^) , 7£ = 0.577215 . . . (4.2) 

e 

Dimensional regularization has many advantages because does not spoil 
the gauge symmetry of QCD and, therefore, simplifies a lot the calculations. 
One could argue that a cut-off procedure would be more physical, since the pa- 
rameter M could be related to some unknown physics at very short distances. 
However, within the QCD framework, both prescriptions are equally meaning- 
less. One just introduces a regularizing parameter, such that the integral is 
well defined and the divergence is recovered in some limit (M ^ 00 or e ^ 0) . 
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Since the momentum-transfer has dimensions, it turns out to be con- 
venient to introduce an arbitrary energy scale /i and write 



^''(4^) r(-e) = V^|^+7iJ-ln47r + ln(V^') + OW} ■ 

Obviously, this expression does not depend on /i; but written in this form one 
has a dimensionless quantity (— inside the logarithm. 

The contribution of the loop diagram in Fig. |l^ can finally be written as 

Kl = ^ab HV"' + 9^'?'^) n(q2), (4.3) 

with 



n{q')^^-Tp |_+^^_ln4^ + ln(-gVM')-- + 0(e) 

Owing to the ultraviolet divergence, this equation does not determine the 
wanted self-energy contribution. Nevertheless, it does show how this effect 
changes with the energy scale. If one could fix the value of Il{q^) at some ref- 
erence momentum transfer gg, the result would be known at any other scale: 



n(g^) = U{qt) - -Tp [^) In {q'/qt) . (4.4) 

We can split the self-energy contribution into a meaningless divergent piece 
and a finite term, which includes the q^ dependence, 

n{q^)^ Mi,{^?) + nR{q^/^?). (4.5) 

This separation is ambiguous, because the finite g^-independent contributions 
can be splitted in many different ways. A given choice defines a scheme: 

■ [7 + - ln(4^) - f ] scheme) , 

An,(Ai') = <( _|£ |i ^2. i (MS scheme), 

-If £ [7 + 7£ - ln(4^)] (MS scheme), 

(4.6) 

If li: In (-^Vm^) (a* scheme), 

nfl(gVM') = {-^i [In (-^Vm') + 1E- ln(4^) - f ] (MS scheme), 
|i [In (-^Vm') - f ] (MS scheme). 
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In the ^ scheme, one uses the value of lU^fi^) to define the divergent part. 
MS and MS stand for minimal subtractionES and modified minimal subtraction 
schemes Ej; in the MS case, one subtracts only the divergent 1/e term, while 
the MS scheme puts also the 7^ — ln(47r) factor into the divergent part. Notice 
that the logarithmic dependence is always the same. 

4-2 Renormalization: QED 




Figure 11: Photon self-energy contribution to e e scattering. 

A Quantum Field Theory is called renormalizable if all ultraviolet diver- 
gences can be reabsorbed through a redefinition of the original fields and cou- 
plings. 

Let us consider the electromagnetic interaction between two electrons. At 
one loop, the QED photon self-energy contribution is just given by Eq. (4.3), 
with the changes Tp — s- 1 and gs — > e. The corresponding scattering amplitude 
takes the form 

^ |i_n(g2) + ...} , (4.7) 

where J'' denotes the electromagnetic fermion current. 

At lowest order, T(g^) ~ a/q^ with a = e^/(47r). The divergent correc- 
tion generated by quantum loops can be reabsorbed into a redefinition of the 
coupling: 

^ {l~Mv,{^?)-Ii^{q'/^?)] ^ {l-nn{q'|^^')] , (4.8) 



afl(Ai ) = ao i 1 + -^M" 



'--'scheme 



"0 - ' (4.9) 



where cq denotes the bare coupling appearing in the QED Lagrangian; this bare 
quantity is, however, not directly observable. Making the redefinition (4.^ 
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the scattering amplitude is finite and gives rise to a definite prediction for 
the cross-section, which can be compared with experiment; thus, one actuaUy 
measures the renormalized coupUng aji. 



The redefinition (4.8) is meaningful, provided that it can be done in a self- 
consistent way: all ultraviolet divergent contributions to all possible scattering 
processes should be eliminated through the same redefinition of the coupling 
(and the fields). The nice thing of gauge theories, such as QED or QCD, is 
that the underlying gauge symmetry guarantees the renormalizability of the 
quantum field theory. 

The renormalized coupling afl;(/i^) depends on the arbitrary scale fi and on 
the chosen renormalization scheme [the constant Cschomo denotes the different 
finite terms in Eq. (^^)] . Quantum loops have introduced a scale dependence in 
a quite subtle way. Both a7j(/i^) and the renormalized self-energy correction 
Tiii{q^ / jj?) depend on /i, but the physical scattering amplitude T{q^) is of 
course ^ independent: (Q^ = —q^) 



4. J^J, ^ \ 1 + + . • • r (4.10) 



T(,^)^-4.J^'J,^^(l + ^ 

1^ Sir 

Q2 \ + 3^ 

The quantity a{Q^) = au{Q^) is called the QED running coupling. The 
usual fine structure constant a — 1/137 is defined through the classical Thom- 
son formula; therefore, it corresponds to a very low scale = — m^. The 
value of a relevant for LEP experiments is not the same: a{M^)-^ — 1/129. 
The scale dependence of a{Q^) is regulated by the so-called (3 function: 



M3^ = a/3(a); /3(a) = /3i- + /Jz f-) + • • • (4.11) 



da 

At the one-loop level, f3 reduces to its first coefficient, which is fixed by (|4.9| ): 

/JP'^^ = I ■ (4.12) 



The differential equation (4.11) can then be easily solved, with the result: 



a{Q^) = ^^^^ . (4.13) 



Since /3i > 0, the QED running coupling increases with the energy scale: 



a{Q^) > q;(Qo) if > Q§; i.e. the electromagnetic charge decreases at large 
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Figure 12: Electromagnetic charge screening in a dipolar medium. 



distances. This can be intuitively understood as the screening due to virtual 
e+e^ pairs generated, through quantum effects, around the electron charge. 
The physical QED vacuum behaves as a polarized dielectric medium. 

Notice that taking = in Eq. (4.10) we have eliminated all depen- 
dences on ln(Q^//z^) to all orders in a. The nmning coupling (4.13) makes a 
resumniation of all leading logarithmic corrections, i.e. 



n=0 



2tt 



In(QVM') 



(4.14) 



These higher-order logarithms correspond to the contributions from an arbi- 
trary number of one-loop self-energy insertions along the intermediate photon 
propagator in Fig. |ll| [1 - UBiq^/n^) + (n^.d^VM^))^ + •■■]■ 



4.3 The QCD Running Coupling 




Figure 13: Feynman diagrams contributing to the renormalization of the strong coupling. 

The renormalization of the QCD coupling proceeds in a similar way. Owing 
to the non-abelian character of SU{3)c, there are additional contributions 
involving gluon self-interactions. From the calculation of the relevant one-loop 
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diagrams, one gets the value of the first /? function coefficient 



S: 



A.|r.«,-lic..^ii^^ (4,15) 

The positive contribution proportional to the number of quark flavours Nf 
is generated by the q-q loops and corresponds to the QED result (except for 
the Tp = i factor). The gluonic self- interactions introduce the additional 
negative contribution proportional to Nq- This second term is responsible 
for the completely different behaviour of QCD: /9i < if Nf < 16. The 
corresponding QCD running coupling, Q!s(Q^), decreases at short distances: 

lim asiQ^) = 0. (4.16) 

Thus, for Nf < 16, QCD has indeed the required property of asymptotic free- 
dom. The gauge self-interactions of the gluons spread out the QCD charge, 
generating an antiscreening effect. This could not happen in QED, because 
photons do not carry electric charge. Only non-abelian gauge theories, where 
the intermediate gauge bosons are self-interacting particles, have this antis- 
creening propertyE2l. 

Although quantum effects have introduced a dependence with the energy, 
we still need a reference scale to decide when a given can be considered 
large or small. An obvious possibility is to choose the scale at which as enters 
into a strong coupling regime {i.e. as ~ 1), where perturbation theory is no 
longer valid. A more precise definition can be obtained from the solution of 



the /3-function differential equation ( 4.11 ). At one loop, one gets 



lnM+^ = InA, (4.17) 

where InA is just an integration constant. Thus, 

"-^^'^= -/3,lnT^VA^) - (4-^') 

In this way, we have traded the dimensionless parameter gs by the dimension- 
ful scale A. The number of QCD free parameters is the same (1 for mass- 
less quarks), but quantum effects have generated an energy scale. Although, 



Eq. ( 4.13 ) gives the impression that the scale dependence of as (/i^ ) involves two 
parameters, and as(/iQ), only the combination (4.17) is actually relevant. 



When /i >> A, as{fi^) 0, so that we recover asymptotic freedom. At 
lower energies the running coupling gets larger; for ^ —> A, as(^i^) —> oo 
and perturbation theory breaks down. The scale A indicates when the strong 
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coupling blows up. Eq. ( |4.18 ) suggests that confinement at low energies is quite 
plausible in QCD; however, it does not provide a proof because perturbation 
theory is no longer valid when — > A. 



4-4 Higher Orders 

Higher orders in perturbation theory are much more important in QCD than 
in QED, because the coupling is much larger (at ordinary energies). The /3 
function is known to four loops; in the MS scheme, the computed higher-order 
coefficients take the valuesE± (Cs = 1.202056903 . . .) 
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(4.19) 



If Nf < 8, /32 < (/Ja < for Nf < 5, while f34 is always negative) which 
further reinforces the asymptotic freedom behaviour. 

The scale dependence of the running coupling at higher orders is given by: 
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5 Perturbative QCD Phenomenology 

5.1 e^e^ Hadrons 




Hadrons 



Figure 14: e+e — >7*,Z* — > hadrons. 

The inclusive production of hadrons in e+e~ annihilation is a good process 
for testing perturbative QCD predictions. The hadronic production occurs 
through the basic mechanism e+e" ^*,Z* — > qq, where the final q-q pair 
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Figure 15: Diagrammatic relation between the total hadronic production cross-section and 
the two-point function n^"'(q). The qq blob contains all possible QCD corrections. The 
dashed vertical line indicates that the blob is cut in all possible ways, so that the left 
and right sides correspond to the production amplitude T and its complex conjugate Tt, 
respectively, for a given intermediate state. 



interacts through the QCD forces; thus, the quarks exchange and emit gluons 
(and q'-q' pairs) in aU possible ways. 

At high energies, where as is small, we can use perturbative techniques to 
predict the different subprocesses: e+e" qq,qqG,qqGG, . . . However, we do 
not have a good understanding of the way quarks and gluons hadronize. Qual- 
itatively, quarks and gluons are created by the q-q current at very short dis- 
tances, X ~ iZ-s/s- Afterwards, they continue radiating additional soft gluons 
with smaller energies. At larger distances, x ~ 1/A, the interaction becomes 
strong and the hadronization process occurs. Since we are lacking a rigorous 
description of confinement, we are unable to provide precise predictions of the 
different exclusive processes, such as e+e^ — s- IGtt. However, we can make a 
quite accurate prediction for the total inclusive production of hadrons: 

(T(e+e^ hadrons) = (7(e+e^ ^ 99 + qqG + qqGG +...). (5.1) 

The details of the final hadronization are irrelevant for the inclusive sum, 
because the probability to hadronize is just one owing to our confinement 
assumption. 

Well below the Z peak, the hadronic production is dominated by the 7- 
exchange contribution. Thus, we can compute the cross-sections of all sub- 
processes e+e^ ^ 7* ^ qq,qqG, ... (at a given order in a^), and make the 
sum. Technically, it is much easier to compute the QCD T-product of two 
electromagnetic currents [J^jj — J^f Qf 9/7^9/]- 

n'^'^iq) ^ijd^x e^«- (0|T {JU^)J':^{Q)^) |0> - {-g>^^q^ + q>^q^) H„n(g') . 
As shown in Fig. the absorptive part of this object (i.e. the imaginary part. 
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which results from cutting — putting on shell — the propagators of the inter- 
mediate exchanged quarks and giuons in all possible ways) just corresponds 
to the sum of the squared moduli of the different production amplitudes. The 
exact relation with the total cross-section is: 



cr(e+e 



hadrons) 



127r Imncmls) 



(5.2) 



Neglecting the small (away from thresholds) quark mass corrections, the 
is given by a perturbative series in powers of as{s): 




F3 + In ( — 



P2 



HI 



^i(^ln(-|+^ln 



+ 0(al) 



(5.3) 



The second expression, shows explicitly how the running coupling as{s) sums 
an infinite number of higher-order logarithmic terms. So far, the calculation 
has been performed to order a^, with the result (in the MS scheme) E3: 



i^i = 1, 



F2 = 1.986- 0.115 iV/, 



F-?. = -6.637- 1.200 iV/- 0.005 A^?- 



1.240 \ ' ' 



(5.4) 



The different charge dependence on the last term is due to the contribution 
from three intermediate giuons (with a separate quark trace attached to each 
electromagnetic current in Fig. |l^). For Nf ~ 5 flavours, one has: 



11 



asis) 



1.411 



12.80 



0(at 



The perturbative uncertainty of this prediction is of order and includes 
ambiguities related to the choice of renormalization scale and scheme. Al- 
though, the total sum of the perturbative series is independent of our renor- 
malization conventions, different choices of scale and/or scheme lead to slightly 
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different numerical predictions for the truncated series. For instance, the per- 
turbative series truncated at a finite order N has an explicit scale dependence 
of order a^'^^. The numerical values of ag and the F„ {n > 2) coefficients 
depend on our choice of scheme (also /3„ for n > 2). 

The prediction for R^+e- (s) above the b-h threshold is compared with 
the data El in Fig. ^, taking into account mass corrections and Z-exchange 
contributions. The two curves correspond to A^~^^ — 60 MeV (lower curve) 
and 250 MeV (upper curve). The rising at large energies is due toJJie tail of 
the Z peak. A global fit to all data between 20 and 65 GeV yieldsEj 

as(35GeV) = 0.146 ±0.030. (5.5) 



5.2 Z —t Hadrons 

The hadronic width of the Z boson can be analyzed in a similar way: 



2 



T(Z^e+e-) ^ — A^-M ^ / '"IM^ 

^ ' I n>l ^ / \ ^/ 

The global factor 

RT = ^'^f[^f'^ (l + few) (5.6) 



contains the underlying electroweak Z — > qfqf decay amplitude. Since 
both vector and axial- vector couplings are present, the QCD correction coeffi- 
cients Fn are slightly different from Fn for n > 2. 

To determine as from Rz, one performs a global analysis of the LEP/SLC 
data, taking properly into account the higher-order electroweak correc.tions. 
The latest as value reported by the LEP Electroweak Working Group E3 is 

as(M|) = 0.119 ±0.003. (5.7) 



5.3 T —>■ Hadrons 

The calculation of QCD corrections to the inclusive hadronic decay of the r 
leptonc3 looks quite similar from a diagrammatic point of view. One puts all 
possible gluon (and qq) corrections to the basic decay diagram in Fig. ^, and 
computes the sum of all relevant partonic subprocesses. As in the e+e~ case, 
the calculation is more efficiently performed through the two-point function 

nf (g) = i J d^xe*^" (0|r(L^(a;)L"(0)^) |0) 

= i^g'-'^q' + cfq") n« (g2) + q^'^ nf {q^) , (5.8) 
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Im{s) 




Figure 16: Integration contour in the complex s-plane used to compute Rt- 

which involves the T-product of two left-handed currents, L'^ = ^"^^{1 — ^^)dg. 
This object can be easily visualized through a diagram analogous to Fig. |l5|, 
where the photon is replaced by a line and one has a tVt pair in the 
external fermionic lines instead of the e+e^ pair. The precise relation with the 
ratio i?T, defined in ( |l.5| ), is: 

The three-body character of the basic decay mechanism, t~ — > Vrudg, 
shows here a crucial difference with e+e^ annihilation. One needs to integrate 
over all possible neutrino energies or, equivalently, over all possible values 
of the total hadronic invariant-mass s. The spectral functions Imn^'^^(s) 
contain the dynamical information on the invariant-mass distribution of the 
final hadrons. The lower integration limit corresponds to the threshold for 
hadronic production, i.e. rriT, (equal to zero for massless quarks). Clearly, this 
lies deep into the non-perturbative region where confinement is crucial. Thus, 
it is very difficult to make a reliable prediction for the s distribution. 

Fortunately, we have precious non-perturbative information on the dynam- 
ical functions n^°'^''(s), which allows us to accurately predict Rr. n^°'^''(s) 
are analytic functions in the complex s-plane except for a cut in the positive 
real axis. The physics we are interested in lies in the singular region, where 
hadrons are produced. We need to know the integral along the physical cut of 
Imn^°'^^(s) = — ^[n^"'^'(s + ie) — n^'^^(s — ie)]. However, we can use Cauchy's 
theorem (close integrals of analytic functions arjE-^ero if there are no singular- 
ities within the integration contour) to expressed Rj- as a contour integral in 
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the complex s-plane running counter-clockwise around the circle |s| = m^: 



Rr = QiTi 4> — ^ ( 1 



n 



(0+1) 



The advantage of this expression is that it requires dynamical information only 
for complex s of order , which is significantly larger than the scale associated 
with non-perturbative effects in QCD. A perturbative calculation of Rr is then 
possible. 

Using the so-called Operator Product Expansion techniques it is possible to 
showtj that non-perturbative contributions are very suppressed [~ (A/ttIt-)^]. 
Thus, Rr is a perfect observable for determining the strong coupling. In fact, 
T decay is probably the lowest energy process from which the running coupling 
constant can be extracted cleanly, without hopeless complications from non- 
perturbative effects. The r mass lies fortuitously in a compromise region where 
the coupling constant is large enough that Rr is very sensitive to its value, 
yet still small enough that the perturbative expansion still converges well. 

The explicit calculation gives E3: 



(5.9) 



where 5ew = 1-0194 aud — 0.0010 are the leading and next-to-leading 
electroweak corrections tZi, and 5p contains the dominant perturbative QCD 
contribution: 



<5p = 
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/5i 
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as (ml) 



+ Oiat) (5.10) 



as{ml) 



5.2023 



asirnr) 



26.366 



+ Oiai). 



The remaining factor (5np includes the small mass corrections and non-perturbative 
contributions. |_. 

Owing to its high sensitivity l3 to the ratio Rr has been a subject 
of intensive study. Many different sources of possible perturbative and non- 
perturbative contributions have beeiuanalyzed in detail. Higher-order logarith- 
mic corrections have been resummedo, leading to very small renormalization- 
scheme dependences. The size of the non-perturbative contributions has been 
experimentally anabized, through a study of the invariant-mass distribution 
of the final hadronsc3; the data implyEj that ^np is smaller than 1%. 
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The present experimental measpxementsS give 5p — 0.200 ± 0.013, which 
corresponds (in the MS scheme) toE3 



a,{ml) = 0.345 ±0.020, 



(5.11) 



significantly larger (lie) than the Z decay measurement (5.7). After evolution 
up to the scale Mz, the strong coupling constant in Eq. (5.11) decreases to 
as{M'^) — 0.1208±0.0025, in excellent agreement with the direct measurement 
at the Z peak and with a similar error bar. The comparison of these two 
determinations of in two extreme energy regimes, mr and Mz, provides 
a beautiful test of the predicted running of the QCD coupling; i.e. a very 
significant experimental verification of asymptotic freedom. 



5.4 



Jets 




Figure 17: Gluon bremsstrahlung corrections to e+e — > qq. 



At lowest order in as, the hadronic production in e+e^ collisions proceeds 
through e~^e~ qq. Thus, at high energies, the final hadronic states are 
predicted to have mainly a two-jet structure, which agrees with the empirical 
observations. At 0{as), the emission of a hard gluon from a quark leg generates 
the e+e~ —>■ qqG transition, leading to 3-jet configurations. For massless 
quarks, the differential distribution of the 3-body final state is given by: 

ctq dxidx2 37r (1 — xi)(l — ^2) ' 

where 

A 2 ^/ 

ao = ^NcY.Q) (5.13) 

^ /=i 

is the lowest-order e+e^ ^ 7* ^ QQ cross-section. The kinematics is defined 
through the invariants s = q^ and = (pi + Pj)'^ = {q — Pk)'^ = s(l — Xk) 
(i,j,k = 1,2,3), where pi, p2 and are the quark, antiquark and gluon 
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momenta, respectively, and q is the total e^e momentum. For given s, there 
are only two independent kinematical variables since 

xi+X2+X3^2. (5.14) 

In the centre-of-mass system, Xi — Ei/Ee = 2_E,;/^/s. 



Eq. (5.12) diverges as xi or X2 tend to 1. This is a very different infinity 
from the ultraviolet ones encountered before in the loop integrals. In the 
present case, the tree amplitude itself is becoming singular in the phase-space 
boundary. The problem originates in the infrared behaviour of the intermediate 
quark propagators: 

xi^l {p2+ Psf = 2 (p2 • P3) ^ ; 

X2^l ^ (Pl+P3)' = 2(pi-P3)^0. 

There are two distinct kinematical configurations leading to infrared diver- 
gences: 

1. Collinear gluon: The 4-momentum of the gluon is parallel to that of 
either the quark or the antiquark. This is also called a mass singularity, 
since the divergence would be absent if either the gluon or the quark had 
a mass (p3||p2 implies S23 = if P2 = P3 = 0). 

2. Soft gluon: pa -> 0. 

In either case, the observed final hadrons will be detected as a 2-jet con- 
figuration, because the qG or qG system cannot be resolved. Owing to the 
finite resolution of any detector, it is not possible (not even in principle) to 
separate those 2-jet events generated by the basic e^e~ — > qq process, from 
e+e^ — > qqG events with a collinear or soft gluon. In order to resolve a 3-jet 
event, the gluon should have an energy and opening angle (with respect to the 
quark or antiquark) larger than the detector resolution. The observable 3-jet 
cross-section will never include the problematic region 0:1^2 ^ 1; thus, it will 
be finite, although its value will depend on the detector resolution and/or the 
precise definition oi jet [i.e. a depends on the chosen integration limits). 

On the other side, the 2-jet configurations will include both e+e^ — > qq 
and e+e^ — )■ qqG with an unobserved gluon. The important question is then 
the infrared behaviour of the sum of both amplitudes. The exchange of virtual 
gluons among the quarks generate an 0{as) correction to the e+e^ ^ qq 
amplitude: 

T[e+e- ^ qq]^To+Ti + --- (5.16) 
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Figure 18: 1-loop gluonic corrections to e+e — > qq. 



where Tq is the lowest-order (tree- level) contribution, Ti the 0{as) correction, 
and so on. The interference of Tq and Ti gives rise to an 0{as) contribution 
to the e+e~ — > qq cross-section. 

We know already that loop diagrams have ultraviolet divergences which 
must be renormalized. In addition, they also have infrared divergences as- 
sociated with collinear and soft configurations of the virtual gluon. One can 
explicitly check that the 0{as) infrared divergence of a{e^e~ qq) exactly 
cancels the one in a{e'^e~ qqG), so that the sum is well-defined: 



a{e e qq) + a{e e 



qqG) 



<7o 



1 + ^ 
TT 



(5.17) 



This is precisely the inclusive result discussed in Sect. 5.1 . This remarkable caxLj 



cellation of infrared divergencesJii actually a general result (Bloch-NordsieckEJ 
and Kinoshita-Lee-NauenbergE3 theorems): for inclusive enough cross-sections 
both the soft and collinear infrared divergences cancel. 




Figure 19: 2-jet configuration. 



While the total hadronic cross-section is unambiguously defined, we need a 
precise definition of jet in order to classify a given event as a 2-, 3-, . . . , or n-jet 
configuration. Such a definition should be free of infrared singularities, and 
insensitive to the details of the non-pcrturbative fragmentation intoJiadrons. 
A popular example of jet definition is the so-called JADE algorithm E3, which 
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makes use of an invariant-mass cut y: 

3 jet <^ Sij = {pi+pjf>ys (Vi,i = 1, 2, 3) . (5.18) 

Clearly, both the theoretical predictions and the experimental measure- 
ments depend on the adopted jet definition. With the JADE algorithm, the 
fraction of 3-jet events is predicted to be: 

where 

Li2(^) ^ ^InC. (5.20) 

The corresponding fraction of 2-jet events is given by i?2 = 1 — i?3. The 
fraction of 2- or 3-jet events obviously depends on the chosen cut y. The 
infrared singularities are manifest in the divergent behaviour of Rz for y ^ 0. 

At higher orders in Ug one needs to define the different multi-jet frac- 
tions. For instance, one can easily generalize the JADE algorithm an classify 
a {j3i,P2, • • • ,Pn} event as a n-jet configuration provided that Sij > ys for all 
i, j = 1, . . . ,n. If a pair of momenta does not satisfy this constraint, they are 
combined into a single momentum and the event is considered as a (n — 1) 
jet configuration (if the constraint is satisfied by all other combinations of 
momenta). The general expression for the fraction of n-jet events takes the 
form: 



Rn{s,y) = 




(5.21) 



with -Rn = 1- A few remarks are in order here: 

• The jet fractions have a high sensitivity to Us [Rn ~ Q^""^]- Although 
the sensitivity increases with n, the number of events decreases with the 
jet multiplicity. 

• Higher-order as(/i^)-' ln'^(s//i^) terms have been summed into q;s(s). How- 
ever, the coefficients cj"''(j/) still contain hi^iy) terms. At low values of 
y, the infrared divergence {y 0) reappears and the perturbative series 
becomes unreliable. For large y, the jet fractions i?„ with n > 3 are 
small. 
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• Experiments measure hadrons rather than partons. Therefore, smce 
these observables are not fuUy inclusive, there is an unavoidable depen- 
dence on the non-perturbative fragmentation into hadrons. This is usu- 
ally modelled through Monte Carlo analyses, and introduces theoretical 
uncertainties which need to be estimated. 

Different jet algorithms and jet variables (jet rates, event shapes, energy 
correlations, . . . ) have been introduced to optimize the perturbative analysis. 
In some cases, a resummation of as{s)"' lii"{y) contributions with m > n has 
been performed to improve the predictions at low y values Ej. 

Many measurements of Og, using different jet variables, have been per- 
formed. All results are in good agreement, providing a consistency test of the 
QCD predictions. An average of measurements at the Z peak, from LET" and 
SLC, using resummed O(a^) fits to a large set of shape variables, givesQ: 

= 0.122 ±0.007. (5.22) 

A recent 0{al) DELPHI analysis of 18 different event-shape observables 
finds an excellent fit to the data, by allowing the renormalization scale to be 
determined from the fit; this procedure gives a more precise value: as{M'^) ~ 
0.117 ±0.003. 

Three-jet events can also be used to test the gluon spin. For a spin-0 gluon, 
the differential distribution is still given by Eq. ( |3.12 ), but changing the x^+x: 
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factor in the numerator to a:|/4. In general, one cannot readily be sure which 
hadronic jet emerges via fragmentation from a quark (or antiquark), and which 
from a gluon. Therefore, one adopts instead a jet ordering, xi > X2 > x^, 
where Xi refers to the most energetic jet, 2 to the next and 3 to the least 
energetic one, which most likely would correspond to the gluon. When X2 — > 
1 {xi 1) the vector-gluon distribution is singular, but the corresponding 
scalar-gluon distribution is not because at that point X3 = (1 — a;i) + (l — 2:2) — > 
0. The measured distribution agrees very well with the QCD predictions with 
a spin-1 gluon; a scalar gluon is clearly excluded. 

The predictions for jet distributions and event shapes are functions of the 
colour group factors Tp = 1/2, Cp = {N^, - l)/(2A^c) and Ca = Nc- These 



quantities, defined in Eq. ( A. 5 ), result from the colour algebra associated with 
the different interaction vertices, and characterize the colour symmetry group. 
If the strong interactions were based on a different gauge group, the resulting 
predictions would differ in the values of these three factors. Since the vertices 
contribute in a different way to different observables, these colour factors can 
be measured by performing a combined fit. The data are in excellent agreement 
with the SU{3) values, and rule out the Abelian model and many classical Lie 
groups. 
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6 Determinations of the Strong Coupling 
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Figure 20: CompilationEfl of as measurements as function of the energy scale. 

In the massless quark limit, QCD has only one free parameter: the strong 
coupling as- Thus, all strong interaction phenomena should be described in 
terms of this single input. The measurements of as at different processes and 
at different mass scales provide then a crucial test of QCD. Obviously, the 
test should be restricted to those processes where perturbative techniques are 
reliable. Moreover, the same definition of as should be taken everywhere; the 
MS scheme is usually adopted as the standard convention. Since the running 
coupling is a function of energy, one can either compare the different deter- 
minations at the different scales where they are measured, checking in this 
way the predicted dependence of the coupling, or use this prediction to 
bring all measurements to a common reference scale where they are compared. 
Nowadays, the Z mass scale is conventionally chosen for such a comparison. 

In order to assess the significance of the test, it is very important to have a 
good understanding of the uncertainties associated with the different measure- 
ments. This is not an easy question, because small non-perturbative effects 
can be present in many observables. In addition, some quantities have been 
computed to a very good perturbative accuracy (next-to-next-to-leading or- 
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Figure 21: Summary of Os measurements, evolved to the scale Mz- 



der), while others are only known at the leading or next-to- leading order; the 
resulting values of as refer then to different perturbative approximations. The 
estimate of theoretical uncertainties is also affected by the plausible asymp- 
totic {i.e. non-convergent) behaviour of the perturbative series in powers of ag. 
Although this is a common problem of Quantum Field Theories, it is probably 
more severe in QCD because the coupling is rather large (at usual energies). 

Fig. ^ summarizes ii the most reliable measurements of the strong cou- 
pling as function of the energy scale. The agreement with the predicted ruiuiing 
of as, indicated by the curves, is indeed very good. Fig. |l] comparesES the 
different measuremenis at the common reference scale Mz. The average af all 
determinations gives 0: 

a^(M|) = 0.119 ±0.002. (6.1) 
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7 Chiral Symmetry 



Up to now, we have only discussed those aspects of QCD which can be analyzed 
in a perturbative way. Thus, we have restricted ourselves to the study of 
scattering processes at large momentum transfers, and inclusive transitions 
which avoid the hadronization problems. The rich variety of strong-interacting 
phenomena governed by the confinement regime of QCD has been completely 
ignored. 

There are certainly many approximate tools to investigate particular as- 
pects of non-perturbative physics; however, rigorous first-principle QCD calcu- 
lations seem unfortunately out of reach for present techniques. Nevertheless, 
we can still investigate some general properties of QCD using symmetry con- 
siderations. 



7.1 Flavour Symmetries 

In order to build the QCD Lagrangian, we made extensive use of the SU (3)c 
colour sym metry, which is the basis of the strong interaction dynamics. The 
Lagrangian ( 3.10| ) has additional global symmetries associated with the quark 
fiavour numbers: 



1- 'Cqcd is invariant under a global phase redefinition of all quark flavours. 

If — ' exp{ie) Qf . (7.1) 

This symmetry is associated with the conservation of the baryon number. 

2. 'Cqcd is also invariant under independent phase redefinitions of the dif- 
ferent quark fiavours. 



Qf — ' expiiOf) Qf . (7.2) 

This symmetry implies the conservation of flavour. 

3. For equal quark masses, there is a larger symmetry under SU{Nf) trans- 
formations in flavour space. 



Qf — > Uff, qj, , U€ SU{Nf) . (7.3) 

This is a good symmetry of the light flavour sector (u, d, s), where 
quark masses can be ignored in flrst approximation. One has then the 
well-known isospin {Nf = 2) and SU (3) symmetries. 
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In the absence of quark masses, the QCD Lagrangian spHts into two 
independent quark sectors, 

1 

4 



£qcd = -7 G^G;:, + z qL^D^qL + i qB^D^qR . (7.4) 



Here, q denotes the flavour (and colour) vector q"^ = {u, d, . . .), and L, R 
stand for the left- and right-handed components of the quarks. Thus, the 
two quark chiralities live in separate flavour spaces (gluon interactions do 
not change the chirality), implying that all previous flavour symmetries 
get duplicated in the two chiral sectors. 



The baryon number symmetry (TA_) is usually called U{l)v, since both 



chiralities transform in the same way. Its chiral replication is the corresponding 
U{1)a transformation: 

qL — > exp(-i6') qL ; qn, — * exp(i6') qp, ■ (7.5) 

This symmetry of the classical (massless) QCD Lagrangian gets broken by 
quantum effects (triangular loops of the type shown in Fig. H, with gluons 



instead of photons); this is the so-called U{1)a anomaly. Although (7.5) is 
not a true symmetry of QCD, it gets broken in a very specific way, which 
leads to important implications. A discussion of the phenomenological role of 
anomalies is beyond the scope of these lectures. However, let me mention that 
this anomaly is deeply related to interesting low-energy phenomena such as 
the understanding of the rj' mass, or the polarized proton structure functions. 

I want to concentrate here in the chiral extension of the old eightfold 
SU{3)v symmetry, i.e. in the global G = SU{3)l SU{3)r symmetry of the 
QCD Lagrangian for massless u, d and s quarks. This larger symmetry is not 
directly seen in the hadronic spectrum. Although hadrons can be nicely classi- 
fied in SU{3)v representations, degenerate multiplets with opposite parity do 
not exist. Moreover, the octet of pseudoscalar mesons [n^K ^rf) happens to be 
much lighter than all other hadronic states. 

There are two different ways in which a symmetry of the Lagrangian can 
be realized. In the usual one ( Wigner-Weyl) , the ground state (the vacuum) 
is also invariant. Then, all physical ^ates can be classified in irreducible rep- 
resentations of the symmetry groupO. Certainly, the hadronic spectrum does 
not look like that in the case of the chiral group. 

There is a second (Nambu-Golstone), more sophisticated, way to real- 
ize a symmetry. In some cases, the vacuum is not symmetric. The hadronic 
spectrum corresponds to energy excitations over the physical vacuum and, 
therefore, will not manifeS|L-.the original symmetry of the Lagrangian. How- 
ever, Goldstone's theoremE3 says that in such a case there should appear a 
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massless scalar for each broken generator of the original symmetry group. If 

the chiral symmetry is realized in this way, there should be eight pseudoscalar 
massless states (Goldstone bosons) in the hadronic spectrum; this is precisely 
the number of states of the lightest hadronic multiplet: the 0~ octet. Thus, 
we can identify the tt, K and 77 with the Goldstone modes of QCD; their small 
masses being generated by the quark mass matrix which explicitly breaks the 
global chiral symmetry of the Lagrangian. 

In the Standard olcctroweak model, the local SU{2)l U{1)y symmetry 
is also realized in the Nambu-Goldstone way. There, the symmetry-breaking 
phenomena is assumed to be related to the existence of some scalar multiplet 
which gets a vacuum expectation value. Since a local symmetry gets (spon- 
taneously) broken in that case, the Goldstone modes combine with the gauge 
bosons giving massive spin-1 states plus the Higgs particle. The QCD case is 
simpler, because it is a global symmetry the one which gets broken. However, 
something should play the role of the electroweak scalar field. Since quarks are 
the only fields carrying flavour, they should be responsible for the symmetry 
breaking. The simplest possibility is the appearance of a quark condensate 



generated by the non-pcrturbativc QCD dynamics. This would produce a 
dynamical breaking of chiral symmetry, keeping at the same time the observed 
SU{3)v symmetry. 



7.2 Effective Chiral Lagrangian 

The Goldstone nature of the pseudoscalar mesons implies strong constraints 
on their interactions, which can be most easily analyzed on the basis of an 
eff'ective Lagrangian. The Goldstone bosons correspond to the zero-energy 
excitations over the quark condensate; their fields can be collected in a 3 x 3 
unitary matrix U{<j)), 



which parametrizes those excitations. A convenient parametrization is given 



V = 



{0\uu\0) = {0\dd\0) = (0|ss|0) < , 



(7.6) 



(7.7) 



by 




(7.8) 
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where 



$(a;) 



A 

71' 



^4 

V2 



V6 



TT 

V2 



V6 



(7.9) 



V 



if- 



2^ , 

Ve/ 



The matrix ?7(0) transforms hnearly under the chiral group, [gL.R G SU{3)L^ii] 



gLQL, QR 



9R QR 



u{4>) 



9RU{cj))gl, (7.10) 



but the induced transformation on the Goldstone fields <j} is highly non-linear. 

Since there is a mass gap separating the pseudoscalar octet from the rest 
of the hadronic spectrum, we can build a low-energy effective field theory 
containing only the Goldstone modes. We should write the most general La- 
grangian involving the matrix U (0) , which is consistent with chiral symmetry. 
Moreover, we can organize the Lagrangian in terms of increasing powers of 
momentum or, equivalently, in terms of an increasing number of derivatives 
(parity conservation requires an even number of derivatives): 



£eff(t/) = I] ^2 



(7.11) 



In the low-energy domain, the terms with a minimum number of derivatives 
will dominate. 

Due to the unitarity of the U matrix, UU'^ = 1, at least two derivatives 
are required to generate a non-trivial interaction. To lowest order, the effective 
chiral Lagrangian is uniquely given by the term 



Co 



■ Tr [a^c/ta^L/] 



(7.12) 



Expanding C/(0) in a power series in $, one obtains the Goldstone's ki- 
netic terms plus a tower of interactions involving an increasing number of 
pseudoscalars. The requirement that the kinetic terms are properly normal- 
ized fixes the global coefficient /^/4 in (7.12). All interactions among the 
Goldstones can then be predicted in terms of the single coupling /: 



C2 



i Tr [9,, $5'"$] 



1 



12/2 



■Tr 



($ df, ^) («> 9^ «>) 



+ ©($7/*). (7.13) 



To compute the tttt scattering amplitude, for instance, is jipw a trivial 
perturbative exercise. One gets the well-known Weinberg result Ea 



T(7r+7r° ^ 7r+^°) = t/f 



(7.14) 
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Similar results can be obtained for tttt An, Gtt, Stt, . . . The non-linearity of 
the effective Lagrangian relates amplitudes with different numbers of Gold- 
stone bosons, allowing for absolute predictions in terms of /. Notice that 
the Goldstone interactions are proportional to their momenta (derivative cou- 
plings). Thus, in the zero- momentum limit, pions become free. In spite of 
confinement, QCD has a weakly-interacting regime at low energies, where a 
perturbative expansion in powers of momenta can be applied. 



It is straightforward to generalize the effective Lagrangian (7.12) to incor- 
porate electromagnetic and semileptonic weak interactions Cj. One learns then 
that / is in fact the pion decay constant / « = 92.4 MeV, measured in 
TT —^ pLV^ decay. The corrections induced by the non-zero quark masses are 
taken into account through the term 

Cm = ^^Tt[M{U + U^)\ , =diag(m„,mrf,m,) , (7.15) 

which breaks chiral symmetry in exactly the same way as the quark mass term 
in the underlying QCD Lagrangian does. Eq. ( |7T^ ) gives rise to a quadratic 



pseudoscalar mass term plus additional interactions proportional to the quark 
masses. Expanding in powers of <I> (and dropping an irrelevant constant), one 
has 



Cm = \v\ {-;^Tr [M^^] + ^Tr [M^^] +0($V/' )} 



(7.16) 



The explicit evaluation of the trace in the quadratic mass term provides 
the relation between the physical meson masses and the quark masses: 

TO^± = (m„ + md)j^ , mlo = (m„ -|- md)jl - e + 0{e^) , 

\v\ \v\ 



where 



\v\ (m„ - rudf 



2/^ (2ms-mu-md) 
Chiral symmetry relates the magnitude of the meson and quark masses to the 



size of the quark condensate. Taking out the common jt'l//^ factor, Eqs. (7.17) 
imply the old Current Algebra mass ratios, 

TO^± m|.+ m^o 3m^ 



md {mu + m,s) {nid + ms) {rUu + nid + ArUs) 
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(7.19) 



and [up to 0{mu — md) corrections] the Gell-Mann-Okubo mass relation 

3m2^ = Ami - ■ (7.20) 

Although chiral symmetry alone cannot fix the absolute values of the quark 
masses, it gives information about quark mass ratios. Neglecting the tiny 0{e) 
effects, one gets the relations 



^ d _ -J<1— — TfL = 12.6 . (7.22) 



2(to„ + rud) 



In (7.21) we have subtracted the pion square mass difference, to take into ac- 
count the electromagnetic contribution to the pseudoscalar meson self-energies; 
in the chiral limit (m„ = = Trig = 0), this contribution is proportional to 
the square of the meson charge and it is the same for _fC+ and tt^. The mass 
formulae (7.21) and ( 7.22| ) imply the quark mass ratios advocated by Weinberg: 



m„ : md : m, = 0.55 : 1 : 20.3 . (7.23) 

Quark mass corrections are therefore dominated by rus, which is large com- 
pared with m„ and m^. Notice that the difference — m-u is not small 
compared with the individual up and down quark masses; in spite of that, 
isospin turns out to be an extremely good symmetry, because isospin-breaking 
effects are governed by the small ratio (m^ — TO„)/ms. 



The interactions in (7.16) introduce mass corrections to the tttt scat- 
tering amplitude ( [7.14 ), 



T(7r+7r" ^ TT+TT-') - (t - mi) / . (7.24) 

Since / ~ /tt is fixed from pion decay, this result is now an absolute prediction 
of chiral symmetry. 

The lowest-order chiral Lagrangian encodes in_a very compact way all the 
Current Algebra results obtained in the sixties cil. The nice feature of the 
chiral approach is its elegant simplicity, which allows to estimate higher-order 
corrections in a systematic way. A detailed summary of the chiral techniques 
and their phenomenological applications can be found in Refs. 42. 

8 Quark Masses 

Owing to confinement, quarks are not asymptotic states of QCD and, there- 
fore, their masses cannot be directly measured. They can only be determined 
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indirectly, through their influenee on hadron properties. Moreover, since they 
are not observable quantities, quark masses need to be properly defined; i.e. 
their values depend on the chosen conventions. 

A possible (intuitive) definilipn is the so-called pole mass: the pole of the 
perturbative quark propagator CJ. However, there is no pole beyond pertur- 
bation theory. Therefore, one expects this quantity to be very sensitive to 
non-perturbative long-distance effects. 

The simplest prescription is to consider in the same way as Og: quark 
masses are just additional couplings in the QCD Lagrangian (3.10). These 
couplings need to be renormalized and, therefore, one gets scale-dependent 
running quark masses. We will adopt the usual MS scheme to define mq(^'^). 

The scale dependence of mq(/i^) is regulated by the so-called 7 function, 



druq 



7i 



— +72 — 

TT V TT / 



(8.1) 



which is known to four loops 



sEll]: 



71 = 2 , 



73 = ttt: 
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•2) 



Here Cn is the Riemann zeta function ( Q2 = 7r^/6, C3 = 1.202056903..., 
(4 = 7r^/90 and Cs = 1.0369277|Sl5 .). The relation between the MS and pole 
masses is known to three loopsBfia. 

Using the (3 function to trade the dependence on /i by as, the solution of 
) is easily found to be 



the differential equation 
niqip^) = niqi^il) exp I 



das l{.as) 



mq{l4) 



Since j{as)/ P{as) is positive, quark masses are smaller at higher energies: 

ruqil GeV2)/m,(Af|) ^ 2.30 ± 0.05 . (8.3) 
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Quark mass ratios are independent of rcnormalization conventions. This 
is the reason why we have been able to fix the hght quark ratios ( [7.23| ) using 
chiral symmetry considerations. Including the next-to-leading 0{^) chiral 
corrections, the ratios of light quark masses have been determined to beCJ: 

TTL 7TL 

— = 0.553 ±0.043; — = 18.9 ± 0.8 . (8.4) 

rrid md 

The quark masses have been estimated with different methods (QCD sum 
rules, lattice, t decays, ...). Although their precise numerical values have 
been rather controversial in the past, some consensus starts to emerge. The 
present status is shown in Table l|. The heavy quark masses are given at the 
quark-mass scale itself, i.e. mq{Tn^), while for light quarks a reference scale 
/io = 1 GeV has been chosen, rflq = mq{l GeV^). 



Table 1: Present values of the quark masses in GeV. mq = rriq (1 GeV^). 



m„ = 0.0046 ±0.0009 


md = 0.0082 ±0.0016 


rus = 0.164 ± 0.033 


■niciml) = 1.2 ±0.1 


mb{ml) =4.2 ±0.1 


mt(TO?) = 165±5 



The |Up and down quark masses have been fixed from the QCD sum rule 
estimate eZI, rf 



rrid 



12.8 ± 2.5 MeV, using the niu/md ratio (8.4). The 



quoted strange quark mass corresponds to 
from the Cabbibo-suppressed t decay width 
agrees with previous QCD sum rules resultst 



ip recent 0{ai) determination 
i msiml) = 119 ± 24 MeV. It 
and jpcentL^ lattice estimates 
(smaller results are found by other lattice groups Ell). The resulting ratio 
fUs/rfid = 20 ± 6 is in nice agreement with the chiral value (8.4). 

The value given for the, charm quark mass reflects the present lattice and 
QCD sum rules estimates t 



The bottom quack mass has been extracted from 
the behaviour of cr(e^e^ — )■ bh) near thresholdc2l, aiid.from a recent unquenched 
lattice calculation of the B meson binding energyEil. When evolved to the Z 
mass scale, it implies mi,{M'^) =pn2.9 ±0.1, in nice agreement with the value 
mb{Ml) = 2.61 ± 0.55 obtainedbl from the 3-jet Z bbG production rateM; 
this provides evidence for the running of quark masses. 

FinaUi^, the top quark mass given in the Table is the value measured at 
FermilabEH, nit = 174.3±5.1 GeV, assuming that it corresponds to a pole mass 
definition. It is in^ood agreement with the result obtained from electroweak 
fits at the Z peakc3. 
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9 Summary 

Strong interactions are characterized by three basic properties: asymptotic 
freedom, confinement and dynamical chiral symmetry breaking. 

Owing to the gluonic self-interactions, the QCD coupling becomes smaller 
at short distances, leading indeed to an asymptotically-free quantum field 
theory. Perturbation theory can then be applied at large momentum transfers. 
The resulting predictions have achieved a remarkable success, explaining a 
wide range of phenomena in terms of a single coupling. The running of as has 
been experimentally tested at different energy scales, confirming the predicted 
QCD behaviour. 

The growing of the running coupling at low energies makes very plausible 

that the QCD dynamics generates the required confinement of quarks and 
gluons into colour-singlet hadronic states. A rigorous proof of this property 
is, however, still lacking. At present, the dynamical details of hadronization 
are completely unknown. 

Non-perturbative tools, such as QCD sum rules and lattice calculations, 
provide indirect evidence that QCD also implies the proper pattern of chiral 
symmetry breaking. The results obtained so far support the existence of a 
non-zero q-q condensate in the QCD vacuum, which dynamically breaks the 
chiral symmetry of the Lagrangian. However, a formal understanding of this 
phenomena has only been achieved in some approximate limits. 

Thus, we have at present an overwhelming experimental and theoretical 
evidence that the SU (3)c gauge theory correctly describes the hadronic world. 
This makes QCD the established theory of the strong interactions. Neverthe- 
less, the non-perturbative nature of its low-energy limit is still challenging our 
theoretical capabilities. 
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Appendix A: SU{N) Algebra 

SU{N) is the group of N x N unitary matrices, UU^ = U^U = 1, with 
det U = 1. The generators of the SU{N) algebra, (a = 1, 2, . . . , iV^ - 1), 
are hermitian, traceless matrices satisfying the commutation relations 

[T'^,r''] = i/'^^'^T", (A.l) 
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the structure constants /"^"^ being real and totally antisymmetric. 

The fundamental representation T° = A°/2 is iV-dimensional. For N = 2, 
A"^ are the usual Pauli matrices, while for TV = 3, they correspond to the eight 
Gell-Mann matrices: 



/010\ 




(o 


-i 0\ 






o\ 




/o l\ 


= 1 




i 










-1 







\0 Oy 




lo 


y 






Oy 




idiocy 






They satisfy the anticommutation relation 



A* = ^ 




{A°,A''} = JjS^^In + 2d"'"=A% 



(A.3) 



where In denotes the iV-dimensional unit matrix and the constants d are 
totally symmetric in the three indices. 

For SU{3), the only non-zero (up to permutations) /"^"^ and d"*"^ constants 

are 



yl56 _ ^^40 _ ^'Jbr _ j. 



(246 



257 



(345 



-f 



367 



V^^ 2' 

^146 ^ ^157 ^ _^247 ^ ^256 ^ ^344 ^ ^355 
^118 ^ ^228 ^ ^338 ^ _2rf448 ^ _2^558 ^ . 



(A.4) 



■2^688 = 



J377 



-2d 



788 



1 



The adjoint representation of the SU{N) group is given by the (A''^ — l)x 



{N^ - 1) matrices (T^) 



be 



if"". The relations 



Tr(rirl) = /-'^/"^'^ = Ca Sab , 
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Cp 



2' 

1 



2A 



Ca = N, 



(A.5) 



define the SU(N) invariants Tp, Cp and Ca- Other useful properties are 

2 

TV 



jahe jcde _|_ jace jdhe ^ jade J^^^ Q . jabe^cde _|_ jace ^dbe _|_ jade ^bce q 



Appendix B: Gauge Fixing and Ghost Fields 

The fields have 4 Lorentz degrees of freedom, while a massless spin-1 gluon 
has 2 physical polarizations only. Although gauge invariance makes the addi- 
tional degrees of freedom irrelevant, they give rise to some technical complica- 
tions when quantizing the gauge fields. 

The canonical momentum associated with G^, n°(a;) = ^£qcd/'5(9oG^) = 
G^^Q, vanishes identically for ii ~ 0. The standard commutation relation 

mx),nUy)]S{x°-y') = ig^-'S^^^x-y), (B.l) 

is then meaningless ior fi = u = 0. In fact, the field G° is just a classical 
quantity, since it commutes with all the other fields. This is not surprising, 
since we know that there are 2 unphysical components of the gluon field, which 
should not be quantized. Thus, we could just impose two gauge conditions, 
such as G" = and VGa = 0, to eliminate the 2 redundant degrees of freedom, 
and proceed working with the physical gluon polarizations only. However, 
this is a (Lorentz) non-covariant procedure, which leads to a very awkward 
formalism. Instead, one can impose a Lorentz-invariant gauge condition, such 
as dfj_G^ = 0. The simplest way to implement this is to add to the Lagrangian 
the gauge-fixing term 

1 

2? 



Cgf = GDid^G^), (B.2) 



where ^ is the so-called gauge parameter. The 4 Lorentz components of the 
canonical momentum 

n^(^) ^ = ^",0--^ 9,0 [d^GD (B.3) 

are then non-zero, and one can develop a covariant quantization formalism. 



Since (B.2) is a quadratic G^ term, it modifies the gluon propagator: 
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Notice, that the propagator is not defined for ^ = oo, i.e. in the absence of the 



gauge-fixing term (B.2). 

In QED, this gauge-fixing procedure is enough for making a consistent 
quantization of the theory. The initial gauge symmetry of the Lagrangian guar- 
antees that the redundant photon polarizations do not generate non-physical 
contributions to the scattering amplitudes, and the final results are indepen- 
dent of the arbitrary gauge parameter ^. In non-abelian gauge theories, like 
QCD, a second problem still remains. 



G 

nnnnmnnrir^ 





< hnHnmrm^ 

q G 

a) b) c) 

Figure 22: Tree-level Feynman diagrams contributing to qg — ► GG. 



Let US consider the scattering process qq — > GG, which proceeds through 
the three Feynman graphs shown in Fig. The scattering amplitude has the 
general form T = J^^ e^^'^ e"^, The probability associated with the scattering 
process, 

V ^ ij-'CJ^-'V 5:4^)4^)* Y.^-^^'* , (B.4) 

A A' 

involves a sum over the final gluon polarizations. One can easily check that 
the physical probability Vt, where only the two transverse gluon polarizations 
are considered in the sum, is different from the covariant quantity Vc, which 
includes a sum over all polarization components: Vc > Vt- In principle, this 
is not a problem because only Vt has physical meaning; we should just sum 
over the physical transverse polarizations to get the right answer. However, 
the problem comes back at higher orders. 

The covariant gluon propagator contains the 4 polarization components; 
therefore higher-order graphs such as the ones in Fig. get unphysical contri- 
butions from the longitudinal and scalar gluon polarizations propagating along 
the internal gluon lines. The absorptive part of these 1-loop graphs {i.e. the 
imaginary part obtained putting on-shell the two gluon lines within the loop) is 
equal to \T{qq — > GG)^. Thus, these loops suffer the same probability problem 
than the tree-level qq — > GG amplitude. The propagation of unphysical gluon 
components implies then a violation of unitarity (the two fake polarizations 
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< — nriTiririn — ^ < 




a') b') c') 

Figure 23: 1-loop diagrams contributing to qq — > qq. 



contribute a positive probability). 

In QED this problem does not appear because the gauge-fixing condition 
d'^A^ = still leaves a residual gauge invariance under transformations satisfy- 
ing D9 = 0. This guarantees that (even after adding the gauge-fixing term) the 



electromagnetic current is conserved, i.e. 



considers the e+e 



8 



5^(eQ*7''^') ^ 0. If one 



77 process, which proceeds through diagrams identical 



to a) and b) in Fig. |2^, current conservation implies fc^J^^ = k'^iJ'^^ = 0, 
where fc^ and k'^, are the momenta of the photons with polarizations A and A', 
respectively (remember that the interacting vertices contained in J^^ are in 
fact the corresponding electromagnetic currents). As a consequence, the con- 
tributions from the scalar and longitudinal photon polarizations vanish and, 
therefore, Vc = Vt- 

The reason why Vc 7^ Vt in QCD stems from the third diagram in Fig. p2| , 
involving a gluon self-interaction. Owing to the non-abelian character of the 
SU{i)c group, the gauge-fixing condition d^C^ ~ does not leave any residual 
invariancefl. Thus, k^,J^'^'' ^ 0. 




d) 




d') 




Figure 24: Feynman diagrams involving the ghosts. 



^ To maintain di_i{Ga)' = after the gauge transformation (3.7), one would need OSda 
gsf°'^'^di_i{59i,)Gc, which is not possible because is a quantum field. 
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Again, the problem could be solved adopting a non-covariant quantization 
where only the physical transverse polarizations propagate; but the resulting 
formalism would be awful and very inconvenient for performing practical calcu- 
lations. A more clever solution consisttS in adding additional unphysical fields, 
the so-called ghosts, with a coupling to the gluons such that exactly cancels 
all unphysical contributions from the scalar and longitudinal gluon polariza- 
tions. Since a positive fake probability has to be cancelled, one needs fields 
obeying the wrong statistics {i.e. of negative norm) and thus giving negative 
probabilities. The magic cancellation is achieved by adding to the Lagrangian 
the Faddeev-Popov termO, 

-Cfp = ^d^^aDf^r , = - gsF'^^'Gi^ , (B.5) 

where (j)"" (a — 1,...,Nq — 1) is a set of anticommuting {i.e. obeying 
the Fermi-Dirac statistics), massless, hermitian, scalar fields. The covariant 
derivative D'^cj)'^ contains the needed coupling to the gluon field. One can 
easily check that diagrams d) and d') in Fig. ^ exactly cancel the unphysical 
contributions from diagrams c) and c') of Figs. ^ a nd p3| , respectively; so 
that finally Vc = "Pt- Notice, that the Lagrangian ( |B.5| ) is necessarily not 
hermitian, because one needs to introduce an explicit violation of unitarity 
to cancel the unphysical probabilities and restore the unitarity of the final 
scattering amplitudes. 

The exact mechanism giving rise to the Cpp term can be understood in a 
simpler way using the more powerful path-integral formalism, which is beyond 
the scope of these lectures. The only point I would like to emphasize here, is 
that the addition of the gauge-fixing and Faddeev-Popov Lagrangians is just 
a mathematical trick, which allows to develop a simple covariant formalism, 
and therefore a set of simple Feynman rules, making easier to perform explicit 
calculations. 
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